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ABSTRACT 


Considering the damage caused to the structures by the earthquake of April 
16, 2016, there is a need to reinforce several structures that were affected. Among 
the types of reinforcement used we can find the use of TADAS type energy sinks 
arranged on inverted V-shaped Steel diagonals. 


In this investigation a proposal for reinforcement for an existing structure 
will be analyzed, this building works as a Hospital that is located in the City of Manta, 
the same that due to the seismic effeets produced had significant damage to its 
structure. 

The use of this type of TADAS type heatsinks placed on inverted V-shaped 
Steel diagonals represents avery efficient reinforcement method which improves the 
structure's behavior since it provides additional stiffness and damping. For the 
analysis of this structure the CEINCI-LAB Computer system will be used applying 
new functions that allow to optimize the data entry time in a simple way. This 
programming is called AHOP3. 

Keywords: New functions of CEINCI-LAB. Spectral Method Reinforced structure 
with “TADAS” heatsinks. 

ANÁLISIS DE ESTRUCTURAS CON DISIPADORES DE 
ENERGÍA TADAS UTILIZANDO CEINCI-LAB 

RESUMEN 

Considerando los daños provocados en las estructuras por el terremoto del 
16 de abril de 2016, se tiene la necesidad de reforzar varias estructuras que fueron 
afectadas. Entre los tipos de reforzamiento utilizados podemos encontrar el uso de 
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disipadores de energía tipo “TADAS” dispuestos sobre diagonales de acero en 
forma de V invertida. 

En esta investigación se analizará una propuesta de reforzamiento para una 
estructura existente, esta edificación funciona como un Hospital que se encuentra 
ubicado en la ciudad de Manta, la misma que debido a los efectos sísmicos 
producidos tuvo daños significativos en su estructura. 

El uso de este tipo de disipadores tipo “TADAS” colocados sobre diagonales 
de acero en forma de V invertida, representa un método de reforzamiento muy 
eficiente el cual mejora el comportamiento de la estructura ya que proporciona una 
rigidez y un amortiguamiento adicional. Para el análisis de esta estructura se 
utilizará el sistema computacional CEINCI-LAB aplicando nuevas funciones que 
permiten optimizar el tiempo de ingreso de datos de una forma sencilla. Esta 
programación se llama AHOP3. 

Palabras Claves: Nuevas funciones de CEINCI-LAB. Método espectral. Estructura 
reforzada con Disipadores tipo “TADAS” 

1. INTRODUCCIÓN 

El objetivo del artículo es presentar nuevas funciones del sistema de 
computación CEINCI-LAB orientadas al reforzamiento sísmico de una estructura de 
hormigón armado de 6 pisos, solo con disipadores de energía por histéresis del 
material, colocados sobre diagonales de acero en forma de “V” invertida. 

Para el análisis sísmico se utilizará un estudio de peligrosidad sísmica 
realizado para a ciudad de Manta para un periodo de retorno de 975 y 2475 años 
aplicando el sismo máximo considerado (MCE). 

El sistema computacional CEINCI-LAB es una herramienta muy útil para el 
análisis de estructuras (Aguiar, 2012; 2014). Con el pasar de los años se ha venido 
implementado nuevas funciones en este sistema computacional logrando así que 
se convierta en un software muy amigable y poderoso para el uso académico como 
profesional. En éste artículo se presentarán las nuevas funciones implementadas 
denominadas AHOP3 al sistema computacional CEINCI-LAB. 

En la Figura 1, se muestra el reforzamiento del parqueadero de la 
Universidad Laica Eloy Alfaro de Manabí, ULEAM, mismo que cuenta con un 
sistema de pórticos de hormigón armado en el cual se ha considerado como método 
de reforzamiento Disipadores tipo “TADAS” sobre diagonales de acero de sección 
circular en forma de V invertida. Aguiar y Del Castillo (2018) 

En el reforzamiento se colocó en la planta baja solo diagonales de acero, 
con esto se logra frenar a la estructura y es factible hacerlo ya que se tienen vigas 
de 40/80 cm; en el segundo piso se colocaron 2 disipadores TADAS en cada uno 
de los vanos reforzados y en los pisos superior 1 TADAS. El reforzamiento se realizó 
en las cuatro esquinas del edificio cómo se muestra en la figura 1. 
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Además de esto se creó un muro de corte en el eje central, con lo que se 
evitó la falla de columna corta, que tenía la estructura, debido a las rampas 
existentes. 



(a) (b) 


Figura 1 Edificio rehabilitado, ubicado en la ciudad de Manta, Parqueadero 
de la Universidad Laica Eloy Alfaro de Manabí ULEAM, a) Vista lateral del 
parqueadero, b) Disipadores TADAS sobre diagonales de acero 

Por otra parte, en la figura 2 a, se muestra como era el Edificio 
Administrativo de la Empacadora Olimar ubicado en la Parroquia Tarqui, de la 
ciudad de Manta, era un edificio de 3 pisos de hormigón armado y una cubierta de 
galvalumen soportada con estructura de acero estructural, como se observa en la 
figura 2 a. Las patologías que produjeron las fallas estructurales son: 
Concentraciones me masas (jacuzzi) en el último piso y columna corta que produjo 
la falla en cabeza de las columnas que se observa en la figura 2 b. 



(a) (b) 


Figura 2 Edificio Administrativo Empacadora Olimar, ubicado en la 
parroquia Tarqui de Manta; a) Edificio antes del terremoto de 2016; b) Colapso del 
“edificio abierto” durante el terremoto del 2016. 

Esta estructura podía ser reforzada, a pesar de tener una gran afectación 
en mamposterías, ya que se comportó de acuerdo a la filosofía de diseño sismo- 
resistente y desempeño estructural adecuado, el cual es: “ante un sismo severo, la 
estructura no debe colapsar y debe ser posible reparar y reforzar a un valor 
razonable”. Esta filosofía de salvar vidas, la mayor parte de los dueños no la 
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comparte y es que no quieren daño por lo que se está cambiando esta forma de 
diseñar estructuras con la de seguridad continua. Sin embargo, por la afectación de 
la columna corta en los pórticos exteriores, la inversión en el mismo sería muy 
costosa y el comportamiento estructural cambiaría debido que se ha desplazado el 
eje inercial de la columna se podría utilizar esta alternativa de reforzamiento y 
disipación de energía para mejorar el desempeño estructural para que no se 
produzcan los mismos daños que se observan en la figura 2b. 

Lo cierto, es que con los programas que se muestran en este artículo, se 
pueden reforzar estructuras (Una alternativa de reforzamiento) y como se verá en 
el siguiente apartado es muy fácil el uso de las nuevas funciones de análisis sísmico. 

2. DESCRIPCIÓN DE LA ESTRUCTURA ANALIZADA 

Se procedió a realizar el análisis estructural de una edificación ubicada en 
la ciudad de Manta, desarrollada en 6 pisos con una altura total de 21.60 metros, la 
cual tiene un sistema de pórticos bidireccionales de hormigón armado con vigas 
peraltadas y losas alivianadas. La estructura presenta una tipología regular tanto en 
planta como en elevación. 

Debido a las afectaciones que tuvo durante el sismo del 16 de abril del 2016 
reforzamiento tuvo la necesidad de realizar un reforzamiento y rigidización de la 
estructura con un sistema de diagonales de acero de sección cuadrada de 240mm 
de lado y 20mm de espesor en forma de “V” invertida conectadas a 2 disipadores 
tipo “TADAS” de 10 placas, localizada en los 2 primeros pisos. Para los pisos 
superiores se utilizarán diagonales de acero de sección cuadrada de 190mm de 
lado y 15mm de espesor en forma de “V” invertida conectada a 2 disipadores tipo 
“TADAS” de 6 placas. El reforzamiento que se indica en este artículo es una 
alternativa didáctica. 

En la figura 3 se indica las secciones de la estructura de hormigón armado 
tanto en los pórticos transversales y longitudinales, además la disposición de las 
diagonales y los disipadores de energía tipo “TADAS”. 



(a) 
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(c) 
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(f) 

Figura 3 Estructura a analizar; a) Vista en Planta; b) Pórticos exteriores 1 y 4 en 
sentido longitudinal; c) Pórticos interiores 2 y 3 en sentido longitudinal, 
d) Pórticos exteriores A y H en sentido transversal e) Pórticos interiores 
B y G. en sentido transversal; f) Pórticos interiores C, D, E, F. en sentido 
transversal. 

En la figura 4, se indica la numeración de los nudos, elementos y grados 
de libertad de la estructura analizada, tanto en el sentido longitudinal como 
transversal, necesarios para el ingreso de datos al programa AHOP3. 



(a) 






































394 


Aguiar F, Palacios L, Palacios B., Olmedo C, Parra Y, Herrera Q. 



(b) 



Figura 4 Pórtico tipo en sentido longitudinal; a) Numeración de Nudos b) 
Numeración Elementos; c) Grados de Libertad 

Para los pórticos en sentido transversal se procede de la misma manera que 
para los pórticos en sentido longitudinal, se enumera los nudos, elementos y 
grados de libertad. 
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3. PROGRAMA AHOP3 

El programa AHOP3 trabaja con las nuevas funciones existentes de 
CEINCI-LAB desarrolladas por Aguiar etal. (2019, 1,2) y utiliza archivos script que 
tienen como finalidad optimizar el tiempo de la introducción de datos que permite el 
cálculo de la rigidez lateral de todos los pórticos existentes para cada sentido, los 
resultados se van grabando en archivos con una identificación única. Estos son 
llamados en dos Scripts adicionales, el primero calcula la matriz total de cada sentido 
y almacena todas las rigideces de todos los pórticos, y la siguiente realiza el cálculo 
de masas, periodos de vibración, factores de participación, desplazamientos 
máximos, derivas, y fuerzas laterales por piso, mediante análisis modal espectral. 
Los Scripts deben ser ejecutados en el orden presentado en la Figura 5. 



Figura 5. Orden de ejecución del programa AHOP3, de acuerdo a la numeración 
indicada. 

A continuación, se describen el funcionamiento y programación de los Scripts. 


3.1. Ingreso de datos y cálculo de KL de cada pórtico 

En la figura 5, se puede observar que primero se ejecuta el programa 
EJ_SolEx.m con el cual únicamente se ingresan los datos de los pórticos en 
sentido X. Después se ejecuta el programa EJ_SolEy.m que sirve para el ingreso 
de los datos de los pórticos en sentido Y. Los datos de entrada que se ingresan 
son para cada uno de los pórticos, en cada uno de los sentidos. 

En el caso de que la configuración estructural indique pórticos iguales en 
los diferentes ejes, la información de entrada se lo realizará una sola vez, como se 
detalla en los siguientes pasos, y se esquematiza en la figura 6: 

a) La entrada de las propiedades geométricas de la estructura 
analizada iniciacon el ingreso de la variable svysp, estas contienen 
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las longitudes de cada uno de los vanos y de las alturas de 
entrepiso, respectivamente. 

b) Seguidamente, se ingresó las propiedades del hormigón y el acero 
que poseen las secciones que se utilizarán en los análisis 
posteriores. 

c) Luego se crea las secciones de hormigón armado mediante el uso 
de la variable SeccionH, esta contiene en la primera columna el 
número del elemento en referencia, la segunda y tercera columna 
posee las dimensiones de base y altura de la sección, la cuarta 
columna contiene el número de elementos que se repetirán y la 
quinta columna contiene un contador de aumento de estos 
elementos, es decir, si se empieza con la numeración en uno y este 
contador es la unidad, se incrementará el número de elementos con 
estas propiedades de uno en uno hasta el número final de 
elementos indicado en la columna cuarta. 

d) En el paso siguiente, se ha desactivado las propiedades de las 
diagonales, debido a que en esta edificación no se ha trabajado con 
estos elementos rigidizadores. 

e) Dado que el programa desarrollado posee la capacidad de generar 
gráficas representativas con las diagonales, es necesario ingresar 
el marco en el cual estos existen, este vector denominado mar debe 
contener siempre un valor, con el fin de no generar errores de 
compilación interna. 

Es necesario realizar los esquemas manuales de numeración de elementos 
y grados de libertad, con el fin de ingresar los vectores elemHor, elemDia y elemTad, 
como se indica en la figura 4. Estos vectores contienen los números de elementos 
de hormigón, acero o diagonales y elementos “TADAS” que conforman la topología 
de la estructura analizada. Es importante señalar que deben numerarse los 
elementos empezando por los de hormigón, seguidamente los elementos de acero 
y finalmente los elementos “TADAS” como consta en la figura 4b. 

Para el caso de esta estructura, se ha propuesto el uso de elementos 
“TADAS” de un tipo en los dos primeros pisos, el resto de los niveles ocupa 
elementos disipadores de diferentes secciones y propiedades; por lo que, es 
necesario ingresar estas magnitudes y almacenarlas en las variables 
correspondientes. 

El programa rigidez, se muestra en el Anexo A y la programación que se 
han utilizado en las rutinas creadas para este caso se los puede obtener en la 
siguiente dirección web: bit.lv/2IJvW0i . 
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Figura 6. Programa EJ_SolEx.m para el pórtico 1 con disipadores de energía 


sobre contravientos. 

Es importante recalcar que el programa rigidez no es fijo para todos los 
casos de análisis, como se indica en la figura 7. Se utiliza la variable elemTad para 
los dos primeros pisos con la modificación del bucle “for” correspondiente que 
permite el cálculo de las propiedades de rigidez equivalente en base a los 
parámetros mecánicos y geométricos de los “TADAS”. Para los siguientes niveles, 
se utiliza la variable elemTadl y el bucle “for” requerido debe ser modificado para 
iniciar en el piso tercero y terminar en el sexto. 
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Figura 7. Modificación de programa rigidez para ingresar diferentes tipos de 
“TADAS” en los pórticos que tienen disipadores de energía sobre contravientos. 

Una vez que se ingresan los datos bajo la secuencia de pasos del 1 al 7 
indicado en la programación de la figura 6, se ejecuta el programa rigidez que 
determina la matriz de rigidez del pórtico con todos los grados de libertad y 
posteriormente se obtiene la matriz de rigidez lateral condensada estáticamente a 
la cual se la ha denominado KL, esta se graba automáticamente en el paso 9. 

El programa rigidez, en forma corta, contiene las siguientes funciones: 

Cálculos iniciales. 

[nv,np,nudt,nudcol,nudvg,nod,nr]=geometria_nudo_viga(sv,sp) 

[X,Y]=glinea_portico2(nv,np,sv,sp,nod,nr) 

[NI,NJ]=gn_portico2(nr, nv, nudt, nudcol, nudvg) 
[GEN]=geometria_nudo_diagonales(nv, np, nudt, mar) 

[NI2,NJ2]=gn_portico(GEN) 

[NI,NJ]=gn_portico3(NI, NJ, NI2,NJ2) 

[CG,ngl]=cg_sismo2(nod,nr,Y) 

[L,seno,coseno]=longitud(X,Y,NI,NJ) 

[VC]=vc(NI,NJ,CG) 

Matriz de Rigidez elementos de hormigón 

[ELEMH]=gelem_portico(SeccionH); 

[KH]=krigidez(ngl, ELEMH, LH, sen H,cosH, VCH, Eh); 
[KA]=krigidez_acero(ngl,ELEMA,LD,senD,cosD,VCD,Es); 

Matriz de Rigidez elementos TADAS 

[diml ,dim2]=size(elemTad); 

Condensación Estática 

KL=kaa- kab* i n v( kbb) * kba 

Como se demostró anteriormente, el programa rigidez no es fijo para todos 
los casos de análisis, por lo que se debe programar de acuerdo con el problema que 
se tenga. Este programa ejecuta las diferentes funciones de CEINCI-LAB y en ellas 
determina la contribución a la matriz de rigidez, por separado, de los elementos de 
hormigón armado; de las diagonales de acero sin disipadores y de los elementos 
con disipadores sobre diagonales. 
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En la figura 8 se indican los datos de entrada para los pórticos 2 y 3 que no cuentan 
con disipadores de energía. 



Figura 8. Parte del programa EJ_SolEx.m p ara e | pórtico 2, que no tiene 
disipadores de energía sobre contravientos. 


Como se habrá podido observar es más fácil la entrada de datos. Se 
recuerda que las instrucciones de la figura 5, vienen a continuación de los datos de 
la figura 4 y todo constituye el archivo EJ_SolEx.m. 

Para los pórticos en sentido Y, se trabaja en forma similar en el archivo 

EJ_SolEy.m 

3.2. Datos para la obtención de la matriz de rigidez en 
coordenadas de piso 

El programa Main1_Matrizrigidez.m es utilizado para sumar la rigidez de 
todos los pórticos en cada sentido, obteniendo la matriz de rigidez total para la 
dirección longitudinal (KLxx) y transversal (KLyy). Posteriormente se ensamblan los 
valores de la rigidez de cada pórtico en una matriz global denominada KE. Ambos 
resultados se graban para ser utilizados en el cálculo modal espectral. 



Figura 9. Programación del programa Main1_Matrizrigidez.m 
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3.3. Análisis modal espectral. 

El programa Main2_AHOP3_eig.m es utilizado para calcular: matriz de 
masas, periodos de vibración, factores de participación, desplazamientos y derivas, 
fuerzas laterales para un espectro de aceleración de un determinado sismo y 
amortiguamiento dado, mediante el análisis modal espectral. La programación se 
presenta en la figura 10. 




Figura 10. Programación Main2_AHOP3_eig.m 


4. RESULTADOS DE LA ESTRUCTURA ANALIZADA 


Una vez ejecutados los programas de ingreso de datos EJ_SOLEx.m, 
J_SOLEy.m y el programa Main1_MatrizRigidez.m, se obtiene las matrices de 
rigidez calculadas tanto para el sentido longitudinal como transversal: 


617166.48 

-299234.73 

40130.05 

-4855.97 

248.82 

48.36 



609271.85 

-294806.86 

39865.48 

-5245.03 

661.90 

174.93 

-299234.73 

314161.25 

-114624.96 

19954.41 

-1061.61 

31.49 
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2.32 
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-79199.99 

7474.99 
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T/m 
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661.90 
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99789.21 

-36176.00 

48.36 

31.49 

-437.76 

4180.59 
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174.93 

2.32 

-651.23 

6755.98 

-36176.00 

29858.56 


Figura 11. Matrices de Rigidez Lateral sentido X, Y. 
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Para realizar el análisis modal espectral, se debe ejecutar el programa 
Main_2_AHOP3eig.m en el cual se obtienen los siguientes resultados: Matriz de 
masa, Periodos de vibración, Aceleraciones obtenidas para el sismo MCE mostrado 
en la figura 12. 



Figura 12. Sismo MCE para Manta. 


Al tratarse de una estructura que cuenta con disipadores de energía se debe 
utilizar un coeficiente de amortiguamiento de 0.27. El análisis sísmico se lo realizará 
mediante el espectro del sismo máximo considerado (MCE), para el cual se ha 
propuesto un factor de reducción a la respuesta sísmica R=4, se determina un 
coeficiente B el cual se obtiene mediante la siguiente ecuación: 



A partir de esta ecuación se determina que el coeficiente B es igual a 1.66; 
para el cálculo de las aceleraciones espectrales, se utiliza los espectros obtenidos 
en un estudio de peligrosidad sísmica para la ciudad de Manta para 975 y 2475 años 
de periodo de retorno para el sismo MCE. 

A continuación, se presentan los resultados obtenidos, tanto para el sentido 
longitudinal como para el sentido transversal: 
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2.6206 

0.1958 
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(b) 

Figura 13. Resultados a) Matriz de Masas b) Periodo fundamental, Aceleraciones 
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En la figura 14 se muestra los valores de las respuestas máximas obtenidas con la 
metodología de combinación CQC para desplazamientos inelásticos, derivas y 
Cortante Basal tanto en el sentido x-x como para el sentido y-y. 

Sentido x-x. 


0.0140 


0.3877 


614.2104 

0.0370 


0.6392 


547.7733 

0.0825 

m Deriva 

1.2657 

% Vt 

478.9601 

0.1247 


1.1714 


394.5771 

0.1696 


1.2474 


291.7427 

°. 1936 


0.6658 


155.9771 


Sentido y-y. 


q ine 


0.0136 
0.0372 
0.0851 
0.1297 
0.1613 
0.1861 
Figura 14. Resu 


Deriva 


0.3767 

0.6565 

1.3315 

1.2371 

0.8782 

0.6906 


Vt 


563.6535 

533.2953 

502.3103 

407.1135 

292.6411 

185.2565 


tados desplazamientos inelásticos, derivas y Cortante Basal tanto 
en el sentido x-x como para el sentido y-y. 


Los valores de deriva son inferiores al dos por ciento (2%) en los dos sentidos, 
siendo la estructura más flexible en el sentido corto de la estructura como se 
presenta en la Figura 15. 



Figura 15. Resultados de derivas para los sentidos X y Y 

La relación H/T para la estructura analizada es 21.6/0.92 = 23.47, 
de acuerdo con Guendelman et al. (2012), establece que para estructuras 
que tengas una relación H/T entre 20 y 40 se consideran como estructuras 
flexibles. 
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5. CONCLUSIONES 

• Como se pudo evidenciar en el análisis de la estructura realizada con 
rutinas que conforman AHOP3 permiten un fácil y rápido ingreso de 
datos que optimiza los tiempos de cálculo para el usuario. 

• Con los resultados obtenidos del análisis se pudo determinar que el uso 
de este tipo de disipadores aporta un amortiguamiento significativo a la 
estructura con el cual se pueden controlar las derivas que se producen. 

• Con las aceleraciones obtenidas del sismo MCE, se pudo determinar 
que las derivas de piso son menores al 2%; tal como se indican en la 
Norma Ecuatoriana de la Construcción NEC 2015 en sus capítulos 
NEC-SE-DS (Peligro Sísmico), NEC-SE-RE (Riesgo Sísmico), 
garantizando que la estructura analizada, de tipología esencial, tenga 
un desempeño de prevención de colapso ante un sismo raro y severo 
con un periodo de retorno de 2500 años y una probabilidad de 
excedencia del 2% en 50 años. (MCE). 
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